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We investigate the influence of morphology and size on the vibrational properties of disordered 
clusters of colloidal particles with attractive interactions. From measurements of displacement 
correlations between particles in each cluster, we extract vibrational properties of the corresponding 
"shadow" glassy cluster, with the same geometric conflguration and interactions as the "source" 
cluster but without damping. Spectral features of the vibrational modes are found to depend 
strongly on the average number of nearest neighbors, NN, but only weakly on the number of 
particles in each glassy cluster. In particular, the median phonon frequency, cjmed, is essentially 
constant for A'^A'^ < 2 and then grows linearly with A'^A*' for A'^A'^ > 2. This behavior parallels 
concurrent observations about local isostatic structures, which are absent in clusters with A^A'^ < 2 
and then grow linearly in number for NN> 2. Thus, cluster vibrational properties appear to be 
strongly connected to cluster mechanical stability (i.e., fraction of locally isostatic regions), and 
the scaling of ojmed with TV A' is reminiscent of the behavior of packings of spheres with repulsive 
interactions at the jamming transition. Simulations of random networks of springs corroborate 
observations and suggest that connections between phonon spectra and nearest neighbor number 
are generic to disordered networks. 

PACS numbers: 63.50.Lm,61.43.Fs,63.22.Kn,64.70.kj,64.70.pv 
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The phase behavior and vibrational properties of en- 
sembles of repulsive particles are determined largely by 
packing fraction [1]. Samples of monodisperse spheres, 
for example, gain structural order and eventually crys- 
tallize with increasing packing fraction [1], giving rise to 
low frequency plane-wave-like phonon modes. In a re- 
lated vein, ensembles of polydisperse spheres gain con- 
tacts with increased packing fraction, leading to vitrifi- 
cation [2] and "soft phonon modes" whose properties de- 
pend on average numbers of interparticle contacts [3] . By 
contrast to these "space-filling" systems, particles with 
strong attractive interactions can form solid-like phases 
at low macroscopic packing fractions [4]. Dilute gels, for 
example, mechanically percolate across large distances 
[5] , and disordered clusters containing relatively few par- 
ticles often self-assemble into structures with large lo- 
cal packing fraction [6]. In this paper we explore how 
cluster morphology and cluster size affect the vibrational 
properties of disordered materials held together by strong 
attractive interactions. New understanding thus gained 
holds potential to elucidate fundamental differences be- 
tween glassy materials composed of particles with attrac- 
tive versus repulsive interactions, to uncover connections 
between vibrational spectra, mechanical stability, and 
the jamming problem, and to discover attributes of a 
disordered cluster that endow it with the properties of 
bulk glasses. 

To date, a diverse collection of disordered systems have 
been observed to display surprising commonality in their 
vibrational properties. Such systems include molecular 
[7], polymer [8], and colloidal glasses [9]. These disor- 



dered solids exhibit an excess of low frequency modes 
that are believed important for their mechanical and 
thermal properties [10]. The low frequency modes also 
appear connected to scaling and mechanical behaviors 
of repulsive spheres near the zero-temperature jamming 
transition. At the jamming point, such disordered pack- 
ings are "isostatic", i.e., they have exactly the number 
of contacts per particle required for mechanical stability; 
if a single contact is removed, the packing is no longer 
stable. Interestingly, marginal stability permits particle 
displacements that maintain isostaticity without energy 
cost; these motions are manifest as low frequency "soft" 
phonon modes [3,11]. When the sample packing fraction 
is increased above the jamming transition, the number 
of contacts per particle increases, the system is stabi- 
lized [12], and the number of soft modes is found to de- 
crease [3]. In fact, the minimum soft mode frequency 
has been predicted to increase linearly with number of 
contacts per particle above the isostatic requirement [3] . 
Recent experiments have found some of these trends in 
thermal packings of repulsive particles [9], but applica- 
tion of such concepts to systems with attractive interac- 
tions has proven difficult. Ensembles of attractive par- 
ticles can achieve isostaticity at arbitrary packing frac- 
tion, and even when they do not have enough contacts 
to be isostatic as a whole, the attractive systems can still 
have local mechanically stable regions [11]. Therefore, 
our study of vibrational properties in clusters of attrac- 
tive particles will provide useful clues about underlying 
mechanisms responsible for the mechanical properties of 
disordered solids. 
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FIG. 1: a. Plot of the temperature-dependent (T= 300.15 
K and 306.45 K) interparticle potential, u(r), as a function of 
particle separation (r) normalized by particle diameter (D). 
The temperature dependence is induced by wetting effects in 
near-critical water-lutidine mixtures, b. The number-fraction 
of locally isostatic configurations per particle, N/ao/N is plot- 
ted versus average number of nearest neighbors, A^'A'^. Solid 
lines are linear fits within two separate regimes, c. Plot sum- 
marizing N and A''A'^ in every cluster, d-f. Experimental 
snapshots of clusters with N= 261, iV]V= 3.95 (d), N= 22, 
NN= 3.91 (e), and N= 22, NN= 4.09 (f). 



Herein we experimentally investigate the influence of 
cluster morphology and size on the vibrational properties 
of disordered clusters of colloidal particles with attrac- 
tive interactions. The disordered clusters with high local 
packing fractions are formed in water-lutidine suspen- 
sions where wetting effects induce fluid mediated attrac- 
tions between micron-sized polystyrene particles. Each 
cluster is characterized by the number of particles it con- 
tains (N) , an average number of nearest neighbors (NN) , 
and a number of local isostatic configurations (N/so). 
Displacement correlation matrix techniques employed in 
recent papers [9] are used to determine the phonon den- 
sity of states of corresponding "shadow" attractive glass 
clusters with the same geometric configuration and inter- 
actions as the "source" experimental colloidal system but 
absent damping [9] . Surprisingly, the spectra and charac- 



ter of vibrational modes depend strongly on NN but only 
weakly on N. The median phonon frequency, uJmed, which 
characterizes the distribution of low and high frequency 
modes, is observed to be essentially constant for NN < 2, 
and then grows linearly with NN for NN > 2. This be- 
havior parallels concurrent observations about local iso- 



static structures, which are absent in clusters with NN 
< 2 and then grow linearly in number for NN > 2. Thus 
cluster vibrational properties appear to be strongly con- 
nected to cluster mechanical stability (i.e., fraction of lo- 
cally isostatic regions), and the scaling oioJmed with NN 
is reminiscent of the behavior of packings of spheres with 
repulsive interactions at the jamming transition. Sim- 
ulations of random networks of springs corroborate ob- 
servations and further suggest that connections between 
phonon spectra and nearest neighbor number are generic 
to disordered networks. 

The experiments employ bidisperse suspensions of 
micron-sized polystyrene particles (Invitrogen) , with di- 
ameters ds = 1.5/xm and cIl = 1.9/^m, and number ra- 



tio 1:2, respectively. Binary mixtures of particles were 
used to minimize crystallization effects. Particles were 
suspended in a mixture of water and 2, 6-lutidine near 
its critical composition, i.e., with lutidine mass frac- 
tion of 0.28. Colloidal particles suspended in this near- 
critical binary mixture, experience temperature depen- 
dent fiuid-mediated repulsive or attractive interactions 
[13]. Interparticle potentials were determined from mea- 
surements of the particle pair correlation function with 
liquid structure theory and image artifact corrections [14] 
(Fig. la). Many different disordered particle clusters are 
created by first suspending particles deep in the repul- 
sive regime (300.15 K), and then increasing the sam- 
ple temperature (to 306. 5K) in situ. Sample tempera- 
ture control was accomplished using an objective heater 
(Bioptechs) connected to the microscope oil immersion 
objective [15]. Particles are confined between two glass 
coverslips (Fisher) with a spacing of '^(1.1 ± 0.05)c?l, 
making the sample system quasi-2D. The glass cell was 
treated with NaOH, so the particle-wall interaction po- 
tential is repulsive at relevant temperatures [16]. The 
global area fraction is '^0.2. Disordered clusters of var- 
ious sizes and shapes self- assemble. Other clusters are 
assembled with aid of laser tweezers [17], either by grab- 
bing particles and adding them to existing clusters, or by 
dragging an optical trap across a cluster and forcing re- 
arrangements. Samples equilibrated for about six hours, 
and video data were collected at a rate of 10 frames per 
second after equilibration for ^ 1000 seconds. 

As noted above, the particle cluster structure is char- 
acterized by several factors, including average number 
of nearest neighbors per particle and number of locally 
isostatic configurations. Neighbors are defined as parti- 
cles spatially separated by less than a cutoff distance, a 
distance equal to the first minimum in the particle pair 
correlation function. Local isostatic regions consist of 
three particles (a, b, and c) that are mutually nearest 
neighbors (i.e., a and b are neighbors, a and c are neigh- 
bors, and b and c are neighbors). Plots summarizing N, 
NN, and N/so for each cluster studied are shown in Fig. 
lb and c, along with experimental snapshots of selected 
clusters (Fig. Id-f). Note that NN tends to increase with 
increasing N for our distribution of cluster sizes, but that 
the increase is not monotonic. The dependence of N/so 
on NN exhibits two regimes. Specifically, Njso/N is for 
NN< 2, becomes non-zero abruptly at NN= 



grows linearly with NN for NN> 2. 
2 as the 



2, and then 
Thus, we identify 



NN= 2 as the "local isostatic" point. 

The vibrational properties of each cluster are extracted 
by measuring displacement correlations of the particles 
within each cluster [9]. Briefly, we extract the time- 
averaged displacement correlation matrix, which, in the 
harmonic approximation, is directly related to the cluster 
stiffness matrix, defined as the matrix of second deriva- 
tives of the effective pair interaction potential with re- 
spect to particle position displacements. These measure- 
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FIG. 2: a. Median frequency, uumed, versus average number of nearest neighbors, A^A''. Two regimes exist. For NN<2, ujmed 
is constant (line is constant fit). For NN>2, ojmed increases linearly with NN (line is a linear fit). In b-d, A''A^ in different 
sample-types is indicated by plotting NN<1.8 as black squares, 2.2<NN<2.8 as red circles, 3.2<A'^A'^<3.8 as blue triangles, 
4:.2<NN as green diamonds. Lines are best constant fits. b. LOmed versus number of particles, N. c. ujrned versus orientational 
order parameter, i/)6. d. ujmed versus average total nearest neighbor spring constant, k. 



ments of the stiffness matrix permit us to construct and 
derive properties of a "shadow" glass system that has the 
same static properties as our colloidal system [9] . Follow- 
ing [18], we expect undamped particles that repel at very 
short-range, and attract due to fluid mediated effects, to 
give rise to solid-like vibrational behavior on time scales 
long compared to particle collision times but short com- 
pared to the time between particle rearrangement events 
[9]. Data were collected over 1000 seconds so that the 
number of degrees of freedom, 8 < 2N < 500, is small 
compared to the number of time frames (> 10 x 2A^) 
[9] . Additionally, we find elements of the stiffness matrix 
to be far above the noise only for adjacent particles, as 
expected. 

Comparing the frequency spectra of clusters with small 
N can be challenging, since not enough modes are present 
for clear identification of a traditional "peak" frequency, 
and since fluctuations can significantly shift mode fre- 
quency. For these reasons we choose to characterize each 
cluster's density of states by its median frequency, Umech 
i.e., we choose the frequency, u>mech such that half of the 
cluster mode frequencies are smaller than cdmed and half 
are larger. Plots of uJmed as a function of average num- 
ber of nearest neighbors, NN, and as a function of total 
number of cluster particles, N (at nearly fixed NN), are 
shown in Fig. 2a and b. 

Surprisingly, LOmed has little correlation with N, ex- 
hibiting linear correlation coefficients of i? = 0.29. This 
effect is even more apparent when the number of near- 
est neighbors is held nearly constant (Fig. 2b). However, 
UJmed depends strongly on the average number of nearest 
neighbors (NN) . We observe two distinct regimes in this 
case. For NN<2, ujmed is constant. For NN>2, oomed 
increases linearly with NN {R = 0.92). Interestingly, 
the dependence of Wmed on NN is very similar to the de- 
pendence of the number-fraction of locally isostatic con- 
figurations per particle, i.e., N/so/N, on NN (Fig. lb). 
These observations suggest that the vibrational proper- 



ties of disordered clusters are strongly dependent on the 
presence of locally rigid elements. Note, we also expect to 
observe a correlation between uJmed and N for our cluster 
distribution; this correlation arises because NN tends to 
increase with N for typical cluster distributions. Thus, 
while we expect the vibrational spectra of a disordered 
attractive cluster to become similar to that of a bulk glass 
as the total number of particles in the cluster increases, 
the underlying mechanism for this effect originates from 
the average number of nearest neighbors in the cluster, 
rather than the total particle number in the cluster. 

Interestingly, the dependence of ujmed on NN is also 
reminiscent of the behavior of hard-spheres in the vicinity 
of the zero-temperature jamming transition [4]. In this 
case, the characteristic frequency, uj* , of excess quasi- 
localized or "soft" modes increases linearly with NN 
when NN>NNc, where NNc is the number of con- 
tacts necessary for isostaticity. Similarly, in our exper- 
iments with attractive particles, cdmed increases linearly 
with NN when NN>2 and when locally rigid elements 
are present. In thermal experiments with repulsive par- 
ticles, UJmed shows a strong linear correlation with uj* 
(R =0.96), and uJmed has a strong linear relationship with 
NN [9]. Thus our observations suggest that similar "jam- 
ming transition" physics may control properties of both 
highly packed "repulsive" glasses, composed of particles 
with repulsive interactions, and disordered clusters com- 
posed of particles with attractive interactions low (over- 
all) packing fraction. 

As a final demonstration of the importance of 
average-number-of-nearest-neighbors versus total- 
number-of-particles in a cluster, consider two clusters 
that look very different (Fig. 1 b and c) but have almost 
the same number of average nearest neighbors {NN). 
These clusters have similar characteristic frequencies 
(i.e., UJmed ^ 3.0 X 10^ and 3.1 x 10^ for Fig. lb and 
c, respectively). On the other hand, two clusters that 
contain the exact same number of particles but have 
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FIG. 3: a. Median frequency, ujmed, plotted versus average 
number of nearest neighbors, NN , from simulations of ran- 
dom matrices, b. Median frequency, Umcd, plotted versus 
number of particles, N, from random matrices with NN— 1 
(closed squares), 2 (open circles), 3 (solid circles), 4 (open tri- 
angles), 5 (solid triangles), 6 (open squares), 7 (solid squares), 
8 (open diamonds), and 9 (solid diamonds). 



different NN (Fig. c and d) possess a set of very different 
characteristic frequencies (aJmed = 3.1 x 10^ and 3.6 x 10^ 
for Fig. Ic and d, respectively). 

As per other calculable cluster properties, uJmed does 
not appear to correlate strongly with many traditional 
structural quantities. The bond orientational order pa- 
rameter, ipQ [19], for example, does not correlate strongly 
with ujjned, whcu NN is held approximately constant 
{R — 0.44) (Figure 2). Similarly, the average stiffness 
between nearest neighbor pairs, k [20], does not correlate 
substantially with Umed, when NN is held approximately 
constant (i? = 0.10). Thus, little correlation exists be- 
tween local order and uJmed- 

Unfortunately, in experimental systems, structural pa- 
rameters cannot be tuned independently. For example, 
clusters with large N and small NN are difficult to find 
or create. Therefore, as a final check on the impor- 
tance of structural quantities other than NN, we ex- 
plored the calculated spectra of randomly generated net- 
works of springs. Random networks of springs, expressed 
as matrices, Kij, were generated following a few sim- 
ple rules that ensure the matrices only contain informa- 
tion about N and NN [21]: (1) Each element, ij, in 
the matrix represents the spring constant between parti- 
cle/coordinate i and particle/coordinate j; (2) the num- 
ber of rows / columns in these symmetric matrices is twice 
the number of particles, representing each coordinate of 
each particle; (3) the number of off diagonal elements 
greater than zero is equivalent to the number of nearest 
neighbors (each particle has the same number of nonzero 
elements); (4) diagonal elements are set such that the 
sum of each row/column is zero, ensuring translational 
invariance. Thus N and NN can be varied independently. 
For every combination of N and NN, 10, 000 random ma- 



trices are generated. uJmed is calculated from the com- 
bination of all generated frequencies (Figure 3). Note, 
many of these networks cannot be duplicated in real ex- 
perimental systems, since nearest neighbor pairings are 
assigned at random and not based on proximity. Never- 
theless, we found that Wmed follows the same trends in 
these simulations as observed in our experiments: uJmed 
has little or no correlation with N (i.e., with NN held 
constant, ujmed changes by less than 5%), but exhibits a 
strong correlation with NN {R > 0.99). Thus uJmed, a 
dominant feature of cluster vibrational spectra, appears 
to be the result of network connectivity, rather than a 
result of specific structure. 

In conclusion, the spectra and character of vibrational 
modes in disordered "attractive" clusters do not depend 
strongly on the number of particles in the cluster, but 
do depend strongly on the average number of nearest 
neighbors and the number of locally isostatic configura- 
tions. The fact that uJmed depends on NN, but not on 
total number and packing fraction, suggests that these 
disordered clusters are a useful model system for network 
glasses (e.g., silica [22]). Network glasses are composed 
of particles (e.g., molecules) that have directional bond 
forming interactions which set A'^A^ [23], leading to the 
formation of solids at low packing fractions. In fact, 
the vibrational [24] and mechanical [25] properties of 
network glasses are believed to depend strongly on A^A^. 
Thus, the disordered clusters composed of particles 
with attractive interactions that we have introduced 
could serve as a convenient model system for network 
glasses and their many applications (e.g., non-crystalline 
semiconductors [26]). 
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